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Abstract 

We use the low-energy effective field theory for holes coupled to the 
staggered magnetization in order to investigate the localization of holes 
on a Skyrmion in a square lattice antiferromagnet. When two holes get 
localized on the same Skyrmion, they form a bound state. The quantum 
numbers of the bound state are determined by the quantization of the col- 
lective modes of the Skyrmion. Remarkably, for p-wave states the quantum 
numbers are the same as those of a hole-pair bound by one-magnon ex- 
change. Two holes localized on a Skyrmion with winding number n = 1 
or 2 may have s- or d-wave symmetry as well. Possible relations with pre- 
formed Cooper pairs of high-temperature superconductors are discussed. 
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1 Introduction 



In the cuprates, high-temperature superconductivity is separated from antiferromag- 
netism by a pseudo-gap regime. It has been conjectured that the relevant low-energy 
degrees of freedom in the pseudo-gap regime are responsible for superconductivity as 
well. Reliably identifying those degrees of freedom by theoretical investigations is a 
highly non-trivial task, because unbiased first principles analytic or numerical calcu- 
lations in microscopic systems such as the Hubbard or t-J model are presently out 
of reach. In lightly doped antiferromagnets, on the other hand, the situation is more 
favorable. First of all, precise numerical simulations of undoped antiferromagnets 0- 
i|, (such as the Heisenberg model) are possible with the loop-cluster algorithm j^], and 
individual doped holes can also be simulated reliably 0, 0]. Second, the low-energy 
dynamics of lightly doped antiferromagnets can be described with a systematic effec- 
tive field theory for magnons and holes. The pure magnon effective field theory has 
been developed in ji-H] and is completely analogous to chiral perturbation theory 



imp 

for the Goldstone pions in QCD 12:]. In the past few years, the systematic effec- 



tive theory for magnons and doped holes has been constructed [13, |14( in complete 
analogy to baryon chiral perturbation theory — the effective theory for pions and 



nucleons |!5l-ll8l|. In contrast to previous attempts to construct effective theories for 



magnons and holes in a square lattice antiferromagnet 19-23), the construction of [14 



is based on a systematic symmetry analysis and provides a complete set of all terms 
contributing to the effective action at leading and sub- leading order. As a result, the 
predictions of the effective theory are exact, order by order in a systematic derivative 
expansion. In particular, the low-energy physics of any lightly doped antiferromagnet 
is described quantitatively once some low-energy parameters (such as the spin stiffness 
or the spinwave velocity of the underlying microscopic system) have been fixed either 
by experiment or by numerical simulations. The effective theory has been used in sys- 
tematic studies of magnon-mediated two-hole bound states [24J and of spiral phases 



25| . Earlier (but somewhat less systematic) studies had been presented in [23|, |26|, [27 



Systematic effective field theories have also been constructed for antiferromagnets on 



a honeycomb lattice [28|, [29J], as well as for lightly electron-doped antiferromagnets 
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Unfortunately, before one enters the high-temperature superconductor or even just 
the pseudo-gap regime, both antiferromagetism and the systematic effective theory 
that describes it break down. While one might expect that one can hence not learn 
anything about high-temperature superconductivity or the pseudo-gap regime from 
the effective theory, the situation may not be entirely hopeless. In particular, the 
effective theory still contains information about what objects may form when the 
theory is about to break down. In this way, we can identify new candidate low-energy 
degrees of freedom for which another effective theory with an extended validity range 
can be constructed. In this paper, we do not yet attempt to construct an effective 
field theory for the pseudo-gap regime. Instead, we concentrate on the identification 
of new low-energy objects that may form when antiferromagnetism is about to break 
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down. 



When antiferromagnetism is weakened, the spin stiffness p s is reduced. In par- 
ticular, if antiferromagnetism is ultimately destroyed in a second order phase transi- 
tion, p s vanishes at the transition. A small value of p s favors topological excitations 
in the staggered magnetization — the order parameter for antiferromagnetism. In 
(2 + 1) dimensions, the topological excitations of the staggered magnetization vector 
are Skyrmions which carry a topologically conserved winding number n G n 2 [S 2 ] = Z 
in the second homotopy group of the order parameter manifold S 2 . The coset space 
S 2 = SU(2) s /U(l)s arises because in an antiferromagnet the SU(2) S spin symmetry is 
spontaneously broken down to the subgroup U(l) s . The possible role of Skyrmions as 
relevant excitations in quantum antiferromagnets has been discussed in several pub- 



lications [3lH45[. Haldane was first to realize that Skyrmions in an antiferromagnet 



are associated with a geometric phase [31]. When Skyrmions proliferate, antiferro- 



magnetic order is destroyed. Read and Sachdev showed that on a square lattice the 
Skyrmion's geometric phase then implies a competing valence bond solid order with 
4- fold degeneracy j32[. The interplay of geometric phases and competing orders has 
been discussed in detail in 43J. The suppression of Skyrmions has been related to 



unconventional deconfined quantum critical points [38l . l39j . It has also been argued 



that a hole localized near a dopant stabilizes a Skyrmion texture in the staggered 
magnetization 33M371] . The analogies between pions in QCD and magnons in ferro- 
and antiferromagnets have been investigated in detail in 40|, |4l|. In particular, it was 
argued that Skyrmions endowed with fermion number 2 may act as preformed Cooper 
pairs of high-temperature superconductivity. Experimental evidence for Skyrmions in 
the lightly doped insulating antiferromagnet La2Cui_ x Li x 04 in an external magnetic 



field has been reported in [44|. Furthermore, the possible role of Skyrmions for the 



superconductivity of Fe based pnictides and chalogenides has been discussed in [45 



In this paper, for the first time, we investigate the localization of holes on a Skyrmion 
using the low-energy effective theory for lightly hole-doped antiferromagnets on a 
square lattice. In particular, we carefully quantize the Skyrmion's collective modes, 
which allows us to unambiguously determine the quantum numbers of single holes as 
well as hole pairs localized on a Skyrmion. 

At the classical level, the mass of a Skyrmion is given by 47ip s . When p s becomes 
small, these excitations hence become energetically favorable. Skyrmions are beyond 
reach of the systematic derivative expansion of the low-energy effective theory for 
magnons and holes. Indeed, when Skyrmions become relevant low-energy degrees of 
freedom, antiferromagnetism as well as the effective theory that describes it are about 
to break down. Still, the effective theory correctly describes the way in which holes 
couple to a Skyrmion excitation in the staggered magnetization order parameter. In 
particular, holes may get localized on a Skyrmion. When two holes get localized on the 
same Skyrmion, they form a bound state which may represent a relevant low-energy 
degree of freedom even when antiferromagnetism gives way to the pseudo-gap phase. 
In particular, such bound states are a potential candidate for preformed pairs whose 
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condensation may ultimately lead to high-temperature superconductivity. In order to 
decide whether this is a viable scenario, in this paper we investigate the symmetry 
properties of Skyrmion-hole bound states in great detail. We find that the p-wave 
states of two holes localized on a Skyrmion with winding number n — 1 transform 
exactly like the two-hole states weakly bound by one-magnon exchange. Two holes 
localized on a Skyrmion may also have s- or d-wave symmetry. Which of these states 
is energetically most favorable depends on the details of the dynamics, and will remain 
a subject for future investigations. 

The rest of the paper is organized as follows. In Section 2 the effective theory 
for the staggered magnetization order parameter is introduced and Skyrmions are 
discussed as classical solutions. The Hopf term is introduced and the collective modes 
of a rotating Skyrmion are then quantized. In Section 3 doped holes are added to 
the effective theory. In Section 4 states of single holes as well as a pair of holes 
(residing in two different hole pockets) localized on a static or rotating Skyrmion 
are constructed and their symmetry properties are investigated. Possible relations to 
the mechanism responsible for high-temperature superconductivity are also discussed. 
Section 5 contains our conclusions. Finally, the case of two holes residing in the same 
hole pocket is investigated in Appendix A. 

2 Skyrmions in the Effective Theory for the Stag- 
gered Magnetization 

In this section we discuss the collective mode quantization of Skyrmions in the low- 
energy effective theory for antiferromagnetic magnons. 

2.1 Effective Action and its Symmetries 

Magnons are the Goldstone bosons of a spontaneously broken spin symmetry SU(2) S 
with an unbroken subgroup U(l) s . Consequently, magnons are described by a 3- 
component unit-vector field e(x) G S 2 in the coset space S 2 = SU(2) S /U(1) S . Here 
x = (xi,X2,t) is a point in (2 + l)-d Euclidean space-time and e(x) represents the 
direction of the local staggered magnetization vector — the order parameter for the 
spontaneously broken spin symmetry. To leading order in a systematic derivative 
expansion, the Euclidean low-energy effective action for the magnons is given by 



Here p s is the spin stiffness and c is the spinwave velocity. The vacuum configuration of 
the effective theory is described by a constant staggered magnetization vector which 
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can be chosen to point in the 3-direction, i.e. e(x) = (0,0,1). Magnons are small 
fluctuations around the vacuum configuration. It should be noted that, in contrast to 
a ferromagnet, antiferromagnetic magnons have a "relativistic" dispersion relation. 

The most important symmetry of the action is the spontaneously broken spin 
symmetry SU(2) S . In the following, global transformations in the unbroken subgroup 
U(l) s will play an important role. Introducing 

e(x) = (sin 6{x) cos (p(x), sin9(x) sin (p(x), cos 0(x)), (2.2) 

these transformations take the form 

T ^e(x) = (sin 6{x) cos(<p{x) + 7), sin 6(x) sm(tp(x) + 7), cos 6(x)). (2.3) 

It should be pointed out that the SU(2) S spin symmetry plays the role of an internal 
symmetry (analogous to chiral symmetry in particle physics). Consequently, its un- 
broken U(l) s subgroup (which is analogous to isospin in particle physics) should also 
be viewed as an internal symmetry. Because of the analogy with isospin, we denote 
transformations in the unbroken subgroup U(l) s by 1(7). 

In addition to the SU(2) S spin symmetry, the effective action has other symmetries 
as well. First of all, due to the relativistic dispersion relation of antiferromagnetic 
magnons, the leading terms in the effective action have an emergent accidental Poin- 
care symmetry which is not present in the underlying Hubbard or t-J model, and which 
will thus be explicitly broken by higher-order terms in the effective action containing a 
larger number of derivatives. The remaining symmetries are the discrete translations 
and rotations of the underlying quadratic lattice. Similar to the spin symmetry, 
the displacements Di by one lattice spacing in the i-direction are also spontaneously 
broken in an antiferromagnet. They act on the staggered magnetization field as 

D *e(x) = -e(x). (2.4) 

Since the shift symmetries Di are spontaneously broken in an antiferromagnet, it is 
convenient to also introduce modified shift symmetries D[ which combine Di with an 
SU{2) S spin rotation g = ia 2 such that 

D *e(x) = ( ei {x),-e 2 (x),e 3 (x)). (2.5) 

Spatial translations by an even number of lattice spacings, on the other hand, remain 
unbroken. Such translations D(xq) by a distance vector x = (:r i, £02, 0) act as 

D ^e{x) = e{x-x ). (2.6) 

Similarly, parametrizing x = (r cos x, r sin x, t), spatial rotations by an angle j3 act as 

°^e(x) = e(0(/3)x), 0(/3)x = (r cos(x + P), r sin(x + P), t), (2.7) 



5 



and a spatial reflection at the xi-axis is represented by 

R e(x) = e{Rx), Rx = {xx,—X2,t) = (r cosx, — t sinx, t). 
Finally, time reversal, which changes the direction of a spin, acts as 

T e(x) = —e(Tx), Tx = (xi,X2,—t) = (r cos x, r sin x, — t)- 
The effective action of eq. fl2.ip is invariant under all these symmetries. 



(2.{ 



(2.9) 



2.2 Classical Skyrmion Solutions 



In particle physics Skyrmions arise as topological excitations in the pion effective field 
theory for the strong interactions 46(, which takes the form of a (3 + l)-d SU(2)l x 
SU(2)r = 0(4) model. In order to distinguish them from their particle physics 
analogs, the topological excitations in the (2 + l)-d 0(3) model are sometimes denoted 
as baby-Skyrmions. For simplicity, here we also refer to them just as Skyrmions. 
Skyrmions are topologically non-trivial classical solutions of the magnon effective 
theory with integer winding number 

n[e\ = J d 2 x Eijt- [diex 9,-e] G Il 2 [S 2 ] = Z, (2.10) 

in the second homotopy group of the sphere 5* 2 . Correspondingly, there is a topological 
current 



^7T 



) • \d v e(x) x d p e(x)} 



(2.11) 



which is conserved, i.e. d^j^x) = 0, irrespective of the classical equations of motion. 
The winding number n[e] = J d 2 x jt(x) is just the integrated topological charge 
density. Under the various symmetries the topological charge density transforms as 



U{l)s 


I[l) Jt{x) : 


= jt(x), 


A 


D 'Jt(x) = 


-Jt(x), 


D[ 


D 'Jt(x) = 


-Jt(x), 


0(/3) 


° W Jt(x) 


= jt(0(P)x), 


R 


R Jt(x) = 


-jt(Rx), 


T 


T Jt(x) = 


-jt(Tx). 



(2.12) 

In particular, the winding number changes sign under the displacements Di and D[ 
as well as under the reflection R and under the time reversal T. 

Let us consider static classical solutions for which the energy 

E[e\ = [ d 2 x ^-die-die (2.13) 
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is minimized. We can write 

< J d 2 x (die±Eijdjexe) 2 

= [ d 2 x (2die ■ die± 2e ii e- (<%e x dje)) = —E[e\ ± 16vm[e], (2.14) 

J Ps 

which implies the Schwarz inequality 

E[e\ >4>irp a \n[e[\. (2.15) 

Skyrmions are minima of the energy in the topological sector with n[e\ = 1, while 
anti-Skyrmions have n[e\ = —1. At the classical level both have a rest energy of 
M.c 2 = 4np s . (Anti-) Skyrmions satisfy the previous inequality as an equality which 
is possible only if they satisfy the (anti-)self-duality equation 

die + (JEijdjex e = 0. (2-16) 

Here a = ±1 distinguishes between Skyrmions and anti-Skyrmions. It is worth men- 
tioning that static (anti-) Skyrmions are mathematically equivalent to (anti-)instan- 
tons of the 2-d 0(3) model 4?|. Using polar coordinates (xi,^) = ^(cosx, sinx), a 
particular (anti-)Skyrmion configuration is given by 

2f n p n 1r n p n u r 2n — p 2n 



e, 



(7,n,p\ 



r >X)=( r 2n + p 2n COS (^), r 2n + M«X)> r 2n + p 2n ) " ( 2 ' 1T ) 



Depending on the sign of a, this configuration describes a Skyrmion or anti-Skyrmion 
of winding number n[e\ = an (with n G M>o) and size p centered at the origin. It 
should be noted that there are many other multi-Skyrmion configurations with differ- 
ent Skyrmions located in different positions. Such configurations would be important 
in investigations of a Skyrmion gas or liquid. Here we concentrate on a Skyrmion 
centered at a single point, possibly with a larger winding number than just n — 1. 
The winding is chosen to arise from the angular ^-dependence which influences the 
rotational symmetry of the Skyrmion and not from the radial r-dependence which 
only influences the finer details of the dynamics. 

The Skyrmion configurations of eq.f l2.17p have a number of zero- modes. In partic- 
ular, their energy remains unchanged when they are shifted to an arbitrary position 
x, when they are spatially rotated by an arbitrary angle /3, or when they are U(l) s 
spin-rotated by an arbitrary angle 7. Interestingly, spatial rotations and U(l) s spin 
rotations act on a Skyrmion in a similar manner, i.e. 

(2f n p n 2r n p n a r 2n — p 2n 

r 2n + p 2n COS (^ + 7), r 2n + p 2n + 7), ^ + p2n 

(2.18) 
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such that 

7(<T7) e-;,„, P (r, x) = 0{l/n) e a ,nA^ x)- (2-19) 

Another zero-mode is related to dilations. Indeed, the energy of a Skyrmion also 
remains invariant under changes of the scale parameter p. A family of Skyrmion 
configurations is obtained by spin-rotating the original Skyrmion of eq. fl2.17p by an 
angle 177 and then shifting it by a distance-vector x such that 

ea, nm (r, x) = D{x) [ /(<T7) e CT ,„, P (r, X )] • (2.20) 

Under the various unbroken symmetry transformations the configuration of eq. (12.201) 
transforms as 



17(1). : 




[r,x) 


= e cr,n,p,2:,7+7o(' r ? X)) 


D[: 


D '*e 


j,x) 


6— <T,n,p,x,'y(j', x)i 


D : 




j,x) 


= Ctr,n,p,x+x ,f( r , X)? 


0((3): 


<: 'CT,n,p,x,') 


j,x) 


= e<r,n,p,0(l3)xa+rip{ r j X) 


R : 


R e 


j,x) 


6— a,n,p,Rx,— yij": X) ■ 



In particle physics Skyrmions play an interesting role in the effective theory for 
the strong interactions. In particular, Skyrmions arise as topological excitations in 
the pion field [if|. While Skyrmions are outside the validity range of the systematic 
low-energy expansion of chiral perturbation theory, they have been used to model 



baryons phenomenologically 48J. Remarkably, the I^S" 3 ] topological winding number 
of the Skyrmions of the strong interactions has the same symmetry properties as the 
baryon number, and is indeed identified with it. The identification of Skyrmions as 
baryons can even be established within the framework of chiral perturbation theory, 
by investigating the electroma gne tic interactions of pions which are affected by a 
Goldstone-Wilczek current 40, 41, 49, 5(|. Since the underlying QCD theory has a 
conserved baryon number current, the conservation of the topological Skyrme current 
is guaranteed beyond the semi-classical regime. 

It is natural to ask whether the winding number n[e\ e IT^S* 2 ] of the Skyrmions in 
an antiferromagnet can also be identified with a conserved quantity of an underlying 
microscopic system, such as the Hubbard model. In particular, in analogy to particle 
physics, one might suspect that the winding number can be identified with the fermion 
number of doped holes. However, this is not the case because the winding number and 
the fermion number have different symmetry properties. In particular, the winding 
number changes sign under a shift D{ by one lattice spacing, while the fermion number 
does not. Hence, unlike in particle physics, in an antiferromagnet the conservation 
of the topological current is not protected by the underlying microscopic dynamics 
and may thus be limited to the semi-classical regime. Interestingly, when holes get 
localized on a Skyrmion, they endow the Skyrmion with their conserved fermion 
number, which may stabilize the Skyrme beyond the semi-classical regime. 



The conservation of the topological current also plays a central role in the scenario 
of deconfined quantum criticality 
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in which dynamically generated gauge fields 
and deconfined spinons are conjectured to appear at a new type of quantum phase 
transition outside the realm of the standard Ginsburg-Landau- Wilson paradigm. In 
fact, the suppression of Skyrmion number violating (so-called monopole) events has 
been argued to change the universality class of the phase transition in the (2 + l)-d 



0(3) model [381 ] . A better understanding of the role of Skyrmions would thus also be 



useful for addressing the issue of deconfined quantum criticality. 



2.3 The Hopf Term 

The integer winding number n[e\ is defined at any instant of time and is conserved for 
topological reasons. Interestingly, there is another topological invariant — the Hopf 
number H[e\ — which characterizes the topology of the order parameter field e(x) as a 
function of both space and time. The integer- valued Hopf number H[e\ G HsfS* 2 ] = Z 
is an element of the third homotopy group of the sphere S 2 . In order to construct the 
Hopf term, it is most convenient to introduce the CP(1) representation 

P(x) = X - (1 + e(x) ■ a) (2.22) 

of the staggered magnetization field. Here a are the Pauli matrices and, as a result, 
P(x) is a Hermitean 2x2 projector matrix that obeys 

P( x )t = p( x ) ; p( x ) 2 = p( x ) ) TrP(x) = 1. (2.23) 

Under a spin rotation g G SU(2) S the matrix P(x) transforms as 

P{x)' = gP{x)g\ (2.24) 

The matrix P(x) can be diagonalized by a unitary transformation u(x) G SU(2), i.e. 

u(x)P(x)u(xy = ±{l + <r 3 )=( 1 Q °Y uu(a:)>0. (2.25) 

We demand that Un(x) is real and positive, which fixes a U(l) s gauge ambiguity and 
uniquely determines u(x) as 



u(x) 



l + e 3 (x) ex(x) — ie2{x) 

y/2(l + e 3 (x)) V ~ei(x) - ie 2 (x) 1 + e 3 (x) 

cos (^9(x)) sin (~6(x)) exp(—iip(x)) 

— sin \^6{x)j exp(iip(x)) cos ^9(x)j 

-jO{x) \ +isin (h 



cos ( -6(x)) +ism ( )-6{x) ) e^x) ■ a, (2.26) 
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where the unit- vector e^x) is given by 

e v (x) = (— sin <p(x), cos ip(x), 0) . (2.27) 
Under a global SU(2) S transformation g, the diagonalizing field u(x) transforms as 

u(x)' = h(x)u(x)g\ u n (x)' > 0, (2.28) 
which implicitly defines the nonlinear symmetry transformation 

h( X ) = exp(ia(,),3) = ( eM f X)) exp( _° a (z)) ) 6 U(l),. (2.29) 

In this way, the global transformations g G SU(2) S of the spontaneously broken non- 
Abelian spin symmetry "disguise" themselves as local transformations h(x) G U(l) s 
of the unbroken subgroup. The global subgroup transformations 1(7) introduced in 
eq.f l2.3p simply lead to a(x) = —7/2. 

The diagonalizing matrix u(x) maps space-time onto the group manifold S 3 of 
SU(2) S . When the (2 + l)-d space-time is also compactified to S 3 , one can relate 
the Hopf number H[e\ G H^S 2 ] = Z to the topological winding number W[u] G 
Tl 3 [SU(2) s } = Tl 3 [S 3 ]=Z, i.e. 

H[e\ = W[u] = 2^Jdt d 2 x e^Tr [(u^d^u) (u ] d v u) {u ] d p u)] . (2.30) 

It should be noted that the evaluation of eq. fl2.30p requires some care. In particular, 
due to the U(l) s gauge fixing Un(x) > 0, u(x) covers only an S 2 subspace of the 
SU(2) S group manifold S 3 . This may seem to imply that the winding number VF[it], 
which counts the number of times the map u(x) covers S 3 , should vanish. However, 
this is not the case because u(x) in eq. fl2.26p is singular at the Skyrmion center where 
^{x) — ~ 1 (i-e. 0(x) = 7r). The singularities which lie on a vortex line encircled 
by e v (x) contribute non-trivially to eq.f )2.30p . Alternatively, one may remove the 
singularities in u(x) by undoing the U(l) s gauge fixing m u (i) > 0, which implies 
that u{x) extends to all of S 3 . Then eq.f l2.30p can be evaluated in a straightforward 
manner. The Hopf term is S'f/(2) s -invariant because 

W[u'} = W[hug ] ] = W[h] + W[u] - W[g] = W[u\. (2.31) 

Here we have used W[g] = 0, which follows because g is constant, and W[h] = 0, which 
follows because the Abelian gauge transformations h(x) G U(l) s are topologically 
trivial in three dimensions, i.e. II3 [C^ (l)s] — n3[5' 1 ] = {0}. 

Under the various relevant symmetries the Hopf number transforms as 



U(l)s 


H[ IM e\ -- 


= m, 




H[ D *e\ = 


H[e\, 


D' t 


H[^e\ = 


H[e\, 


Oifi) 


H[°^e\ 


= H[e\, 


R 


H[ R e\ = 


-m, 


T 


H?e\ = 


-m- 



(2.32) 
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The Hopf term gives rise to an additional factor exp(iOH[e\) in the Euclidean path 
integral with G being the anyon statistics angle. In systems with reflection or time- 
reversal symmetry, the value of is hence limited to or tt. As we will see, in these 
cases Skyrmions are quantized as bosons or fermions, respectively. In systems without 
reflection and time-reversal symmetry, arbitrary values of G are allowed, and then the 
Skyrmions may have any (neither integer nor half-integer) spin. By investigating field 
configurations in which two Skyrmions interchange their positions, one can also show 
that Skyrmions pick up a phase exp(i6) and thus obey anyon statistics 15 lL I t should 
be noted that the Hopf term is expected to be absent in doped cuprates |3jJ,l32j, 52l - 
while it is known to be present, for example, in quantum Hall ferromagnets |40|, 
59(. In order to keep the discussion as general as possible, we will include the Hopf 
term, although in the cuprates one expects 9 = 0. 



2.4 Collective Mode Quantization of the Skyrmion 

Let us now consider the collective mode quantization of the Skyrmion. The main goal 
is to understand the quantum numbers of the quantized Skyrmion, first of all in an 
undoped system. It should be pointed out that Skyrmions in an undoped antifer- 
romagnet are heavy objects whose pair-creation is suppressed at low temperatures. 
When antiferromagnetism is weakened by hole doping, the Skyrmion mass is reduced 
and, in addition, the holes may lower their mass by getting localized on a Skyrmion. 
This favors Skyrmion formation in doped antiferromagnets. A central goal of this 
paper is to understand the quantum numbers of the Skyrmion-hole bound states. In 
this subsection, we consider the collective mode quantization of a Skyrmion in the 
undoped system. 

In order to perform the collective mode quantization, we consider the zero-mode 
parameters pit), x(t), and j(t) as functions of time. We now evaluate the Euclidean 
action (including the Hopf term) for a time-dependent Skyrmion and (after a some- 
what lengthy but straightforward calculation) we obtain 

S[e a , n ,pw\ + i&H[e a>ntPtX „] = J dt (Owe 2 + ^-x 2 + ^j^p 2 + 7 2 + in ^^ ■ 

(2.33) 

Here the Skyrmion's rest energy is given by 

/^ n 2 2n ? ,2n-2 
^ (r 2n + p 2n)2 = 4 ^ ( 2 - 34 ) 

which confirms that for self-dual solutions the Schwarz inequality of eq. (I2.15P is obeyed 
as an equality. For n > 1 the Skyrmion's inertia against dilations takes the form 

DM= ^/<f> 2 ;::^ . (2.35) 
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For n — 1 the integral is logarithmically infrared divergent. In a finite volume or in a 
system with a finite density of Skyrmions, the infrared divergence may be regularized 
because the volume available to each Skyrmion becomes effectively finite. Indeed such 
effects are known to arise in the instanton gas of the 2-d 0(3) model [i^, 56]. In this 
paper, we do not attempt to decide whether the same happens in the (2 + l)-d 0(3) 
model that is relevant here. We just regularize T>(p) by an infra-red cut-off R which 
may or may not be infinite such that for n = 1 



V{p) 



4irp s 



R 



dr- 



2r 3 



M log 



R 2 + P 2 



R 2 



[r 2 + p 2 ) 2 

Finally, the moment of inertia of the Skyrmion is given by 

Ar 2n p 2n V(p)p 2 



R 2 + p 2 



d x- 



r 2n _|_ p 



2n\2 



n- 



(2.36) 



(2.37) 



which is affected by the same infrared divergence as T>(p). Hence, although the 
Skyrmion has a finite mass (and can thus undergo translational motion), in the limit 
of an infinite infra-red cut-off R it has an infinite moment of inertia X(p) and can thus 
not rotate. 



From the Euclidean action of eq. 02. 33!) we read off the real-time Lagrange function 



as 



M. 2m 

2 2 



•2 , P 2 -2 
P + —1 



e 

— n — 7 

2vr [ 



Mc 2 . 



(2.38) 



In the next step, we consider the canonically conjugate momenta 

V{p)p 2 i 



9L ... dL . . . dL 

Pi = W =M*i, Pp =gr=V(p)p, Pj =- 



(2.39) 



It should be noted that, at the classical level, the 0-term is suppressed because relative 
to p 7 it is of order H (which we have put to 1). The canonically conjugate momenta 
lead to the classical Hamilton function 



H = Piii + p p p + p 7 7 - L = Mc 2 + —j- + 



2M ' 2V(p) 
The momentum p i: the spin p 7 , and the energy 



2 n 2 ( 

v ^T 2 V^ n Y, 



(2.40) 



E 



2V(p) 



2 n 2 ( 
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n 



2V(p)p 2 







(2.41) 



of the coupled rotational and dilational motion are conserved quantities. The last 
equality determines the size p(t) of the Skyrmion as a function of time 



P(t) 



dp 



P(0) 



2E 



n 



V(p) V(p) 2 p 2 



p 7 + n 



e_ 

2^ 



-1/2 



(2.42) 
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When the Skyrmion is rotating (i.e. when p y + ^ 0), centrifugal forces lead to an 
unlimited increase of p(t). 

Upon canonical quantization the momentum pi and the spin p 7 turn into the 
operators 

Pi = -id Xi , p 7 = -id 7 , (2.43) 

while the classical Hamilton function % turns into the quantum mechanical Hamilto- 
nian 

H = Mc 2 - J-dl - —L= ( dl + -d ] -J= - (dL + in^-] 2 . 



(2.44) 

The collective mode wave function of a Skyrmion or anti-Skyrmion with winding 
number an, momentum pi, and spin p 7 = am G Z takes the form 

y p ,a,n,m(x, P, l) = expfexi) eicp(iamj)ijj(p). (2.45) 
The dilational part of the wave function solves the Schrodinger equation 



1 /„, 1„\ 1 n 2 ( e^ 2 



di + -d , == + — : = m + na— 



4>{p)=Ei,{p) ) 



(2.46) 

which may again lead to an instability of a rotating Skyrmion against unlimited 
increase of its size p. As we will see later, localized holes prevent the increase of p 
and thus stabilize the Skyrmion. 

In the presence of the Hopf term the spin operator of the Skyrmion (which is 
analogous to isospin in particle physics) is given by 

I = *(p 7 + n^j = a ^-zo> 7 + n^j . (2.47) 

The state ^ p>(Tjn>m {x, p, 7) hence has the "isospin" 

I^p,a,n,m(x, P, l) = {m + an ^\ ^ p,(T,n,m(x , P, l) ■ (2A8) 

In particular, for = the "isospin" is an integer, while for G = tt it is a half-integer 
for odd n. 

Let us also investigate the quantum numbers of the Skyrmion with respect to 
spatial rotations. As a consequence of eq. (12. 191) . the angular momentum J is given by 

J = anl = n ( p 1 + n— J = n ( — id 7 + n—) , (2.49) 

V 27 7 V 27 7 

J^p,cr,n,m(x, P, l) = n (am + %p,a,n,m(x, P, l)- (2-50) 
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such that 



Hence, for = the Skyrmion has integer angular momentum and thus is a boson, 
while for G = 7r the angular momentum is a half-integer and the Skyrmion is a 
fermion. Interestingly, in (2 + 1) dimensions it is possible to have particles of any 
(neither integer nor half-integer) angular momentum — the anyons which arise for 
6 ^ or 7T. 

By construction, the Skyrmion state is also an eigenstate of the momentum op- 
erator with eigenvalue j?;. Under the modified shift symmetries D[ and under the 
reflection R the Skyrmion state transforms as 

UR^p^mix, P, 7) = ®Rp,-*,n,m(%, P, l), (2.51) 

where Rp = (pi, — P2) is the spatially reflected momentum. Here Ud>. and Ur are 
unitary transformations representing the corresponding discrete symmetries in the 
Hilbert space of the collective modes of the Skyrmion. It should be noted that shifted 
or reflected Skyrmions (which have a = 1) are actually anti-Skyrmions (with a = — 1). 



3 Effective Action for Doped Holes 

In order to make the paper self-contained, in this section we review the main features 
of the effective field theory constructed in [14| which couples doped holes to the 
staggered magnetization order parameter. 



3.1 Nonlinear Realization of the SU(2) S Symmetry 

In order to couple holes to the staggered magnetization order parameter, a nonlinear 



realization of the spontaneously broken SU(2) S symmetry has been constructed in [13 
The global SU(2) S symmetry then manifests itself as a local U(l) s symmetry in the 
unbroken subgroup. This is analogous to baryon chiral perturbation theory in which 
the spontaneously broken SU(2)l ® SU(2)r chiral symmetry of QCD is implemented 
on the nucleon fields as a local SU(2)l = r transformation in the unbroken isospin 
subgroup. 

The definition of the nonlinear realization of the SU(2) S symmetry is based on the 
diagonalizing matrix u(x) defined in eq. fl2.26p . which transforms as 



u(x) = u(x)*, (3.1) 



under the modified displacement symmetry D[. Introducing the traceless anti-Hermi- 
tean field 

v^(x) = u(x)d ll u(x)\ (3.2) 
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one obtains the following transformation rules 



SU(2) 8 : 








A': 




(s) 




: 




(s) 


= EijVjiOx) 


R : 


% 


(*) 


= vi(Rx), 



°v t {x) = v t (Ox), 

R v 2 ( x ) = -v 2 {Rx), R v t (x) = vt(Rx). (3.3) 



Writing 



v^{x) = iv~{x)cr a 



v ti x ) = v l( x ) Tiv 2 Jx), 



(3.4) 



the field v^(x) decomposes into an Abelian "gauge" field v (x) and two "charged" 
vector fields v^(x). 



Using eq. fl2.20p . for a Skyrmion e atnt p t0 ^(r, X ) centered at x = one obtains 

■in 



vl(r,x) 
tf(r, x) 
vf(r, X) 

vf{r,x) 



onp 



r (j,2n _|_ p2n^ 
2n 

— cosx, 



onp 



r(r 2n + p 2n ) 

J2n 



op 



7, 



y>2n __|_ pin 

^ \2n+ p 2n eX P(^ \{n + 1) X + j}), 
r 2n + p 2n [(/l + l) X + 7]), 

or n p n 

^2n _|_ p2n 



exp(=Fi(x(?2x + 7))7- 



(3.5) 



In principle, when holes get localized on a Skyrmion, they affect the radial profile of 
the Skyrmion. Here we neglect this effect and concentrate on symmetry considerations 
which are independent of such details of the dynamics. 



3.2 Hole Fields and their Transformation Properties 



As discussed in detail in |14J the holes are described by Grassman- valued fields ip±(x). 
Here / G {a, (3} is a flavor index which specifies the momentum space pocket in 
which the hole resides, and the subscript ± denotes the spin of the hole relative to the 
direction of the local staggered magnetization. Under the various relevant symmetries 
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Figure 1: Elliptically shaped hole pockets centered at (±^,±7^). Two half-pockets 
combine to form the pockets for the flavors f = a, (3. 



of the underlying antiferromagnet on a square lattice, the hole fields transform as 



SU{2) S : 




x)' 


= exp(±ia(x))ip±(x) , 


U{l) Q : 




(x) 


= exp^i^tp^x), 


D 'i ■ 




(x) 


= ± exp(ik{ a)^£(x), 


: 




(x) 


= T^i(Ox), °i>l(x 


R : 




(x) 


= ^i(Rx), R il>l{x 



^(Rx). 



(3.6) 



The U(1)q symmetry is just fermion number, while k a = and k 13 = (^, — 

(with a being the lattice spacing) point to the centers of the two hole pockets illus- 
trated in figure 1. It is interesting that in the effective theory momentum indices of 
the underlying microscopic dynamics turn into internal flavor quantum numbers. 



3.3 Effective Action for Holes coupled to the Staggered Mag- 
netization 

Based on the above symmetry properties, the leading and sub-leading terms of the 
effective action for an antiferromagnet on a square lattice have been constructed sys- 
tematically in [lij]. Here we restrict ourselves to the leading terms. We also make 
the simplifying (but somewhat unrealistic) assumption that the momentum-space hole 
pockets have a circular shape, which enables us to perform large parts of the following 
calculations analytically. It would be straightforward to take into account the more 
realistic elliptic shape of the hole pockets, but this would require some numerical 
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work. Here we concentrate foremost on the symmetry properties of holes localized 
on a Skyrmion on which the simplifying assumption of spherical hole pockets has no 
effect. The total action of the coupled system including doped holes then takes the 
form 

S[^MA = J d 2 x dt {d^ ■ d t e+ j 2 d t e- d t ej 

(3.7) 

Here M and M' are the rest energy and the kinetic mass of a hole, and A is the 
hole-one- magnon coupling constant. The sign aj is + for / = a and — for f = fl. 
The covariant derivatives are given by 

D^ f ± (x)= %±ivl(x)]il)i(x). (3.8) 

Remarkably, the Shraiman-Siggia term in the action, which is proportional to A, 
contains just a single (uncontracted) spatial derivative. Due to the nontrivial rotation 
properties of flavor, this term is still 90 degrees rotation invariant. Due to the small 
number of derivatives it contains, this term dominates the low-energy dynamics. In 



particular, it alone is responsible for one-magnon exchange between hole pairs [14J, [24 
as well as for potential spiral phases in the staggered magnetization order parameter 
251 ] . It is interesting to note that a similar term is absent in lightly electron-doped 



antiferromagnets 30j, such that spiral phases do not arise in these systems. 



4 Hole Localization on a Skyrmion 

In this section we apply the effective theory of the previous section to the localization 
of holes on a Skyrmion. First, we consider the localization of a single hole first on 
a static and then on a rotating Skyrmion. Then the localization of two holes on the 
same Skyrmion is considered, and the symmetry properties of the resulting two-hole 
bound states are analyzed. 



4.1 Single Hole Localized on a Static Skyrmion 

As we have seen, the moment of inertia X(p) of a Skyrmion with n = 1 is loga- 
rithmically divergent in the infra-red. Unless the divergence is regularized due to a 
finite spatial volume or the presence of other Skyrmions, the Skyrmion then cannot 
rotate. In the interest of analytic solubility, and because we want to focus on sym- 
metry aspects, we will no longer consider the translational and dilational motion of 
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the Skyrmion. Instead, we fix the Skyrmion center at the origin x = and we fix 
the Skyrmion size to a constant p. As we will see later, in the presence of holes, the 
energy of the Skyrmion-hole bound states is minimized for a particular value of p. 



The wave function of a single hole localized on a Skyrmion takes the form 



H,n(r, x) 



(4.1; 



Omitting the constant rest energy M of the holes, which just amounts to a constant 
energy shift, the corresponding Hamiltonian resulting from the action of eq. (13.7j) is 
given by 



Hi 



h{ + h{_ 

H f _ + HL ] ' 



-[d t + tvf(x) 12 



2M 



2M> 



dl + -d r + -[d Y + i 



onp 



2n \ 2 



2 \ \ r^2n _|_ pin 



A(v+(x) + a f vt(x)) 
V2Aoof — — — exp I — to 



yln _|_ pin 



[n + l)x + 7 + Of 



7T 



Hi 



A(Vj (x) + <TfV 2 (x)) 

^ Aaa f r 2n + p 2n eX V\ /a 

1 -[d t -ivKx)\ 2 



{n + l)x + l + °f 



71' 

4J 



2M 



2M' 



dl + -d r + — [d 



onp 



2n \ 2 



f r \ ^ v^ n -\- p^ n 



■ (4-2) 



Using the explicit form of vf(x) and vf(x) for the Skyrmion of eq.f l3.5p and making 
the ansatz 



\[// 

<r,m+,m- 



_(r)exp (io-[m + x- | -o-/f]) 
.(r) exp (zcr[m_x + | + 



(4.3) 



,m_ 

with m_ — m + = n+1, after some algebra one obtains the radial Schrodinger equation 



H r 1p m +,m- 

with 



Hr++ H r+ - 
H r — 1_ i7 r J \ 1pm + ,m-,-( r ) 



m+,ni— ,+ 1 



E 



m+ ,m_ V / m+ 



ip m+ ,m.(r), (4.4) 



Hr++ 

H r +- 



1 



2AP 

i?5 — 4- - 



m + + 



np 



2/i 



^>2n _|_ p2n 



^2n _|_ p2n 



1 



2M' 



c> 2 + -d r - — m_ - 



np 



2)i. 



j~2n _|_ p2n 



(4.5) 
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It should be noted that the resulting radial Schrodinger equation is the same for 
Skyrmions and anti-Skyrmions as well as for both flavors / = a, (3. Interestingly, for 
odd n and m_ = — m + — (n + l)/2, the two equations decouple. The equation that 
leads to a localized hole takes the form 



1 



2M' 



dt + -d r - - 



1 fn + 1 




v / 2Anr n - 1 p n 



„2n 



+ P- 



2n 



where ip(r) is the linear combination 



(4.6) 
(4.7) 



For even winding number n, on the other hand, the two equations do not decouple. 
In the following we will be most interested in Skyrmions (or anti-Skyrmions) with 
winding number n — 1. 

In this paper, we concentrate on the symmetry properties of holes localized on a 
Skyrmion, not paying much attention to finer details of the dynamics. Hence, here 
we do not solve the radial equation, which would be straightforward using numerical 
methods. Still, we want to obtain at least a rough estimate for the ground state 
energy of a hole localized on a Skyrmion. For n — 1 the radial Schrodinger equation 
takes the form 



d z r + -d r ) + V{r) 



ip{r-) = Eip(r), 



2M' 

with the potential given by 

V(r) , , 

1 ; 2M'(r 2 + p 2 ) 2 r 2 + p 2 

At short distances, the potential can be approximated by a harmonic oscillator 



- V2A- 



(4.8) 



(4.9) 



^approx^) 




r 2 + C(r 4 ), 



(4.10) 



and hence, in a rather crude harmonic approximation, the ground state energy takes 
the form 



En 



V2A 
P 



1 



m ,2 p 4 1 My 



+ 



M'A 2 x (^Jx 2 + 2V2x - \/2j , 



1 



x 



M'Ap 
(4.11) 

Minimizing the energy as a function of x yields x 3 + 3\/2x 2 + Ax = y/2, which is solved 
by 



x = 



+ 



1/3 



+ 



+ 



-1/3 



V2^ 0.271 



1 



0.271M' A' 



(4.12) 
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Figure 2: The potential V(r) (solid curve) together with its harmonic approximation 
(dashed curve) and the corresponding ground state energy (dotted line). 



This shows that the presence of the hole explicitly breaks the scale invariance that 
led to the dilational instability of the pure Skyrmion. The resulting bound state with 
the strongest binding energy has 

E = M'A 2 x (^\/x 2 + 2V2x - \/2j « -0.135M'A 2 . (4.13) 

The potential V(r) is shown in figure 2 together with its harmonic approximation and 
the corresponding ground state energy E . The figure implies that the true ground 
state energy is smaller than the harmonic approximation suggests. 



4.2 Single Hole Localized on a Rotating Skyrmion 

In this subsection we consider a single hole localized on a rotating Skyrmion. When 
the moment of inertia X(p) diverges (as it is the case for n — 1 and R = oo) the 
fixed orientation 7 of the Skyrmion explicitly breaks the U(l) s symmetry and the 
analysis of Section 4.1 applies. Here we assume that X(p) = V(p)p 2 /n 2 is finite. This 
is actually the case when the feedback of the localized hole on the radial structure 
of the Skyrmion is taken into account. When X(p) is finite, the Skyrmion can rotate 
and thus 7 becomes a dynamical variable. The 7-dependent terms in the Lagrange 
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function for the rotational motion are given by 



2n 2 



(4.14) 



f=a,P 
»=+,- 



Using eq. fl3.5p . the momentum canonically conjugate to 7 thus takes the form 



p 7 



n 



— n- 

2 2n 



d x 



(7- 



P 



2n 



V f 



=a,/9 



(4.15) 



which leads to the corresponding Hamiltonian 

1 



with the Berry gauge field 

A,= I d*x £ 



»=+■■ 



OR g\fyf _ ji^- 



(4.16) 



(4.17) 



Combining the results, one sees that while the off-diagonal elements of the Hamil- 
tonian ( 14. 2 p remain the same, the diagonal elements receive additional contributions 
such that now 
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(4.18) 



21 



We now make the ansatz 



\Sjf 

a,m+,m-,m 



1pty,m+ ,m— ,m,+ ( 



ia[m 



exp (ia [m + x — °7§]) ex P 
r) exp (ia [m_x + c/f ] ) exp(i<r( 



m 



m J 

(4.19) 

with m_ — m + — n + 1. In order to ensure 27r-periodicity of the wave function in the 
variable 7, m must now be one half of some odd integer. This is in contrast to the 
rotating Skyrmion without a hole that was discussed in Subsection 2.4, for which m 
was an integer. The radial Schrodinger equation is then given by 



m+,m-,m\ 



r++ 



r+- 



i>a. 



m+ ,m- ,m 



_(r) 



^a,m+,m-,m^a,m+,m-,m(j')' (4.20) 

In this case, the four matrix elements of the radial Hamiltonian H r take the form 
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(4.21) 



4.3 Symmetry Properties of a Single Hole Localized on a 
Skyrmion 

Let us again consider the spin operator (which generates an internal symmetry and 
is thus analogous to isospin in particle physics) 

f -iad 1 + an% + \ \ , . 

which commutes with the Hamiltonian, i.e. [H* , I] = 0. The wave function ^,m+,m_,m 
is indeed an eigenstate of /, i.e. 

I H,m + , m -,Jr, X , 7)=(m + an^j H,m + , m .,Jr, X , l)- (4.23) 
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Since m is half of an odd integer, the rotating Skyrmion with one hole localized on it 
has half- integer spin (or "isospin"), at least for vanishing anyon statistics parameter 
9 = 0. 

The various symmetries such as the displacements D[ and D' 2 , the 90 degrees 
rotation O, as well as the reflection R, act on the wave function 

Kn(r, x ,i) = ( *r + S r '* ,7 0> (4 - 24) 

of a single hole localized on a rotating (anti-) Skyrmion with winding number an as 
follows 



Di H,n(r, X, 7) = exp(^fa) ( *' 7) 

ff ' nl ' X ' 7J V *L-(r,X + f,7-n|) 



2 ' 

/ 



yf (rv 7 )=( *W(»-.-X,-7) 



(4.25) 



For energy eigenstates this then implies 

D **a,m+,m_,m( r )X)7) = ^/ exp(ifcfa)tf{. ff _ m _ _ m+ _ m (r, x, 7) , 
0$ ",m + , m .,m( r .X,T) = exp (ia[m + + m_ - 2 - 2nm]^) #J m+ , m _ )m (r, X, 7), 

°^, m+ , m „, m ( r ,X,7) = - exp ^[m + +m_-2nm]^) *£m +1 m_,m( r > X> 7), 
fl \l/ a fr y 7) = ^ y 'y! 

ff,m^,m_,mv > A) // ^ — Gr,m-|_,m_,mv ) A) (/) 
^ a,rn+,m-,rn{ r i X> 7) a-,m-)-,m_,m('"> X) 7) ■ (4.26) 

It should be noted that for ^ or 7r, the reflection symmetry i? is explicitly 
broken by the Hopf term. Assuming appropriate phase conventions for the radial 
wave functions, in the considerations of the shift symmetries D^, we have used 

rri— ,— m,+ y) . .m .m.— v) ; 1p—m-,—m+,—m, — VVre+,m_,m,+ (^)5 (4.27) 

which follows from the behavior of eq.f l4.2ip under the replacement of m + — > m' + = 
— m_, m_ — > m'_ = —m + , and m — > m! = —m. It is worth noting that after this 
replacement the constraint 

m'_ — m' + = —m + + m_ = n + 1 (4.28) 

remains satisfied. 
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4.4 Schrodinger Equation for a Pair of Holes of Different Fla- 
vor Localized on a Rotating Skyrmion 



Let us now consider bound states of two holes localized on the same Skyrmion. Both 
a hole of flavor a and another hole of flavor /3 can occupy the same single-particle 
ground state in a Skyrmion. For holes of the same flavor this would be forbidden 
by the Pauli principle. Since we are most interested in the lowest energy states, we 
consider two holes of different flavor. The case of two holes with the same flavor is 
discussed in Appendix A. The Hamiltonian for two holes of different flavor a and f3 
is given by 

E = E a + E? + E\ (4.29) 

where E a and E 13 are the Hamiltonians for a hole of flavor a and (3, respectively. 
Explicitly one has 
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(4.30) 
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= H V T( X ) - v7 ( x )), 



-(d t -ivf(x)) 2 , Hi 



n 



2V{p)p 2 
n 2 

mp)p 2 

n 2 
n 2 

'2V(p)p 2 



) , ■ 6 
a 7 + in- ta 

Z7l 



Jin 



2n 



r z a n + p 

e . P 2n 

o-y + in- ta 

2n 



cL + in- — h ia— — 

2tt r - 1 - n 



r 2 a n + p 1 

2n 

r^ n + p 

gL + in h ia— — 

7 2tt r 2n + p 2n 



— l<7 



ia 



P 



2n 



rf + p 



In 



Jin 



rj n + p 



,2n 



— ia- 



P 



In 



e 



r 2 g n + p 

Jin 



,2n 



+ ia- 



rf + p- 



,2n 



(4.31) 
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We now make the following ansatz for a two-hole energy eigenstate 



a/3 



o,m°^_ ,m" ,m^_ ,vnT_ ,m 



( ip 

— (Tip 



a a p (9 . (r a , r« I exp ta 



m\Xa + m p + xp 



exp(ia(m — 1) 



\ 



^m^m^mLm^-^a,^) eX P [™ 



u,m°^_ ,m^_ ,m^ + ,ml ,m 



<7,m™ ,m™ ,m' + ,m'_ ,m, 



m%Xa + m'_Xf3 



m a _x a + m p + xp + f 



exp(iam / y) 
J exp(iani'y) 

__(r a , r p ) exp (V m^Xa + m p _xp ) exp(m(m + 1)7) 



(4.32) 



Again, this solves the Schrodinger equation only if mL — — n + 1. As for the 
Skyrmion without holes, in this case, m is again an integer. The resulting radial 
Schrodinger equation then takes the form 



ff.mx ,m_ ,vnr. ,mr ,m 



^ ^0 m a yyjP 

with 



E 



r °7 



m °t ,m™ ,m 



/ "0 




s 


s 






\ 




,m« 






,m,+-( r "' 








,m« 




S 










,m« 






,m, ( rQ? ' 




/ 



The radial Hamiltonian is given by 



with 



( H« 



r++ 



TTa 
rL r- + 





m = 
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sl r-+ 





V 






H r+- 
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rra 





TTa 
I2 r+- 







Tja 

ii r 







Tja 





H?__ 
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H r+ _ 
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Hf__ 






















H r -+ 


H r . 


/ 



H 1 

r 



T + + + + 







TTl 

S1 r+— 









tt-t 

S1 r-+— 









H2 



\ 
) 



(4.34) 



(4.35) 



(4.36) 
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The matrix elements of the fermionic part of the radial Hamiltonian are 

2" 



H. 



f 



1 



-O r _| 



2M' 



1 



2n 



d rf H d r , 5- mi + q~~ T - 

s r-f f r j \ r j + P 



„n-l 



nr, p" 
#/_+ = \/2 A „ f - , 

r .2n + 2n ' 



2M' 



^ + -^ ,2 



1 / * np 
i „p _ 



2n 



r jn _|_ p2r. 



(4.37) 



while the rotational Skyrmion contributions are given by 



H 1 



TT-y 



n 



H2 



2V(p)p 2 
n 2 

n 2 
2V(p)p 2 

n 2 



m + era- 1 — 



2ra 



n 2n 



27T r^ + p 2 " r 2n + p 2n 



p 2n p 
m + an- — h 



2n 



m + an h 



'2n rl n + p 2n r 2n + p 2n ) ' 

Q n 2n n 2n 



2vr rl n + p 2n r 2n + p 2n ' 



„2n 



9 , 1 , l"" P 

m + an- — h 1 H — — h 



2n 



27T ' ' r2™ + p 2 ™ r 2n + p 2n 



(4.38) 



4.5 Symmetry Properties of a Pair of Holes with Different 
Flavors Localized on a Skyrmion 

It is worth noticing that the spin operator /, which commutes with the two-hole 
Hamiltonian H, is given by 

/ -iad^ + an^- + 1 \ 

* -iad y + an§; 

-iad^ + an^ 

\ -iad 1 + an^ - 1 J 

(4.39) 



such that 



m + an^-)ty alB p p (r a ,Xa,rp,Xfai)- ( 4 - 40 ) 

27J- I a,m",m'±,m l ^ r ,m_,m^ K H v ' 

Since m is an integer, as expected, for 6 = the state with two holes localized on a 
Skyrmion has integer spin (which plays the role of "isospin" ) . 
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The symmetries D[, O, and R act on a general two- hole wave function 



( * Xa,rp,Xp,l) \ 

Ki,- + ( r a,Xa,Tp,Xp,l) 

V (r a ,Xa,rp,Xp,l) ) 



(4.41) 



as follows 



^",n( r a, Xa, Tp, X/3, 7) = exp(i(fc? + fcf )a) 



-R,T r a 



Xa,Tp,Xp,l) \ 

K,n-+(r a ,Xa,rp,xp,l) 
v-{r a ,x a ,rp,xp,l) 

°,n,++( r a)Xa) r /3)X/3)7) / 



/ ~ fo, X/3 + f , ?"a, Xa + f , 7 ~ ™f ) \ 



_+(r^, X/3 + f , r a , X« + f , 7 - ™f ) 
^af n ,+-(rp,xp + ^,r a ,Xa + f ,7 - nf) 

V Cl— ( r /3>X/3 + f,^a,Xa + f,7-«f) / 

/ -Xp,r a , -Xa, -7) \ 

Cn,-+( r /3, -X/3, r a , -Xa, -7) 

-X/3, r a , -Xa, -7) 

V ,—irp, ~Xp, r a , -Xa, -7) / 



(4.42) 



It is straightforward to show that for the two-hole energy eigenstates this implies 



D'^aP 



<j,m^_ ,m_ ,m + ,m_ ,m 

D '^ a " a a P P 

a,m^_ ,m_ ,m^_ ,m_ ,m 

0, Tr aP 



p (r a ,Xa,rp,xp,7) = y 



aP 



<j,—m°i,—m°l_,—m'~_,—m'j_,—m 



(r a ,Xa,rp,Xp,l) = -® a t 



K^_, m 'i, m tJ r ^ *«' r ^ 7) = exp ( m[m + + m - 



ir a ,Xa,rp,Xp,l), 

_l_ , m 

Xa,rp,xp,l), 
mn + 1]— J 



R-q,ap 



a a P 

cr,rrv± ,m_ ,m , ,m_ ,m 



X Knl,m-,rnP + , m t,J r °> X <*> ^ X ^' 
,aP 



T a , Xa, Tp, X p, 7) = *Z,rnP + , m i, m l, m «_J r °> X ™ X ?> l] (443) 



Here we have assumed an appropriate phase convention for the radial wave function 
ib a a p p (r„,r«). In the context of the shift symmetries D'- we have used 

t a.m. ,m a ,m, ,m ,m> a, p j j % 



/3 


^Sr*,rp) 


= lb 


/3 /3 
+ ,"»- 








, m , + Sra,rp) 


= lb 
"r-a 






= lb 



-m" ,— ,— rrt^ ,— ,— m,++ 

-m_ ,— m\_,— m,-\ — v l J ' ' 

_ m /3,_ m ^,_ m ,_ + (ra,r /3 ), 



%°L, -m£,- 
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These relations follow from the symmetries of the radial Schrodinger equation f)4.33p . 
Similarly, in the context of the rotation O we have used 

^o-,m^,m^,m^,m^,m,++( r / 3 ' r °) = ^ o,m? + ,m^_ ,m°[_,m™ ,m,++ ( rQ ' T ' 
^o-,m« ,m^,m^,m^,m,-+( r ^' r «) = ^ o,m? + ,m p _ ,m,+- ( r ° ' T P) ' 

a Q /3 /3 . (r8,r a )—lb 13 13 a a i(V/»,7\}), 

r tT,m^,m_,mij_,m_,m,H — V M' "/ r <r,m\_ ,m_ ,m^_,m_ ,m. — h v ' rV' 

^Ar.m" ,m° ,m, v .3' ^V*) ^a,m,^ ,mP_ ,m™ ,m° ,m, vc*) ^/?) " (4.45) 

Finally, in the context of the reflection symmetry R we have used 

a a a p , , (rg,rVv) = "0 a a a . . (/v^rs), 

^<x,m™ ,m" ,m, — h v ' <r,m^ ,m" ,m" ,m.H — ' '"V') ' 

^<r,m™ ,m™ ,»n,H — v ' ^a) ^—cr,m^ ,mP_ ,m" ,m" ,m. — (- ' ^ ft) ' 

■0 a a /? (r8,r a ) = ib /3 i3 a a (r a ,r 3 ). (4.46) 

The relations in eq. (14.461) follow from the symmetries of the radial Schrodinger equa- 
tion (I4.33P for = 0. For G ^ or 7r, the Hopf term explicitly breaks the reflection 
symmetry. 



4.6 Comparison with Two-Hole States Bound by One-Mag- 
non Exchange 



In [14| states of two holes bound by one-magnon exchange in a square lattice antifer- 
romagnet have been investigated in great detail. Here we summarize as well as extend 
some of the relevant results. In the rest frame, the Schrodinger equation for two holes 
of flavor a and (3 takes the form 



M 1 



E l ^ • (4.47) 



The components *&i(r) and ^(r) are probability amplitudes for the spin-flavor com- 
binations and respectively. The potential 

27rp s r z 

couples the two channels because magnon exchange is accompanied by a spin-flip. 
Here r = f + — r_ is the distance vector between the two holes of spin + and — and (p 
is the angle between r and the x-axis. Magnon exchange is attractive between holes of 
opposite spin, and hence magnon-mediated two-hole bound states are invariant under 
the unbroken subgroup U(l) s . We make the ansatz 

*i(f) ± * 2 (f) = R(r) X ±(<p). (4.49) 
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Figure 3: Angular wave functions ceo(<p, M'A 2 /47rp s ) (solid curve) and 
sei((p, M' A 2 / Anp s ) (dashed curve) as well as the angle- dependence cos(2<p) of the 
potential (dotted curve) for a pair of holes with flavors a and (3 (M'A 2 /4np s = 1.25 ). 



For the angular part of the wave function this implies 

" "3^ ± ^003(2^^) = -X X± W). (4.50) 
This is a Mathieu equation whose solution with the lowest eigenvalue — Ai is given by 



X±0) = ^ ce ° yf^ 



M'A 2 



4vrp s 



M'A 
4vrp s 



2\ 2 



+ (9(A S 



(4.51) 



The first excited state and its eigenvalue — A2 is given by 



M'A 2 
47rp s 

I / 7T M'A 2 
se 1 [ip- 



2 4vrp s 

M'A 2 1 /M'A 2X 2 

A 2 = -l + - + - 

47rp s 



7T 



cei <p 



M'A 2 
4vrp, 

C(A 8 ). 



1 

64 



M'A 2 
4vrp s 



(4.52) 

8 V 47rp s J 64 V 4vrp s / v ' y ' 

For small A, A2 < 0, which (as we will see) implies that the corresponding two-hole 
state is unbound. For M 1 A 2 /Aixp s > 0.908046, on the other hand, Ai, A 2 > 0, such that 
then both states are bound. The periodic Mathieu functions ce (<p, M' A 2 / An p s ) and 
sei(<p, M'A 2 /Anp s ) j6o| are shown in figure 3. The corresponding radial Schrodinger 
equation is given by 

~d 2 RAr 



dr 2 



r dr 



^Ri(r) = M'EiRi(r), i G {1, 2}. (4.53) 
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Figure 4: Probability distribution for two holes with flavors a and (3. Left panel: the 
ground state with p-wave symmetry. Right panel: excited states with s- or d-wave 
symmetry, but with identical probability densities (M'A 2 /47rp s = 1.25, tq = a). 

The short- distance repulsion between two holes can be incorporated by a hard core of 
radius r , i.e. we require Ri(r ) = 0. The radial Schrodinger equation for the bound 
states is solved by a Bessel function 

R i (r)=A i K v ( y /M'\E ik \r), k = 1,2,3,..., v = % (4.54) 
The energy (determined from K^i^J M^E^Tq) = 0) is then given by 

E ik ~ -{M'rlY 1 exp(-2irk/y/Xi) (4.55) 

for large n. Magnon exchange mediates weak attractive forces that lead to a small 
binding energy. 

The two lowest energy states with angular part x+i^p) an d X-if) are degenerate 
in energy. Linearly combining the two states to two eigenstates of the rotation O, one 
obtains 

The corresponding probability density is illustrated in figure 4 (left panel). While the 
probability density seems to resemble d x 2_ y 2 symmetry, unlike for an actual d-wave, 
the wave function is suppressed, but not equal to zero, along the lattice diagonals. In 
fact, as one operates on the states \?±(f) with the 90 degrees rotation O, one obtains 
the eigenvalues ±z, which shows that they actually have p-wave symmetry. 
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Under the discrete symmetries D[, O, and R, the ground states \I/±(r), which are 
bound by magnon exchange, transform as 



D' , T ,1 



i?i(r) 

i*i(r) 
i2i(r) 



X^O) ±*X-0) 
X\{<P) T iX-{<p) 

X+if) TiX-(v) 

xl(v + f)±<xi(v+f) 

xi(^) ±*x+(v) 
-xi(v) ±*x+M 

x+O) ±*x-H 

x+O) T«X-0) 



±i*i(r), 



(4.57) 



It should be noted that in [14| there are two typos in the last line of the previous 
equation for the reflection symmetry R (eq.(6.20) in [Hj]). 

Remarkably, the magnon- mediated two- hole ground states ^±(r) transform ex- 
actly as the two- hole states localized on a rotating Skyrmion with n — 1, provided 
that we associate *&±(r) with the corresponding two-hole-Skyrmion wave function 



m 



m_ 



aiXai r p-,XPil) with the quantum numbers cr = ±, m + 
1, and m = 0. Indeed, according to eq. (14.431) one obtains 



m , 



£ ^*£-l,l,-l,l,o( r «> X«> r /? ? X/3, 7) 



±.-l,l,-l,l,o( r «>Xa)^>X^>7) 



±.-l,l,-l,l,o( r «>Xa)^>X^>7) 



**±?-i,i,-i,i,o( r «> Xa, r/?, 7) 



-l,l,o( r «, Xa, X/3, 7), 
-*^-l,l ) -l,l,o( r «> Xa, Xf3, 7), 

,-i,i,o( r «) Xa, r p , xp, 7), 
-i,i,o( r «) Xa, rp, xp, 7)- (4.58) 



Just as the magnon-mediated bound states, these states are also invariant under U (l) s 
and they have fermion number 2. One may argue that the two-hole-Skyrmion states, 
in addition, have Skyrmion number as a conserved topological quantum number. 
However, as we discussed before, Skyrmion number has no analog in the underlying 
microscopic Hubbard or t-J models and is just an accidental symmetry of the effective 
theory. We thus conclude that, in their ground state, two holes bound by magnon 
exchange indeed have exactly the same quantum numbers as two holes localized on a 
rotating Skyrmion with n — 1. This implies that these sets of states may evolve into 
each other upon doping. In this way, two holes weakly bound by magnon exchange at 
small doping may evolve into a strongly correlated preformed pair of holes localized 
on a Skyrmion. However, as we have just seen, theses bound states actually have 
p-wave symmetry. 
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Let us also consider the excited states, bound by magnon exchange, with angular 
part x+iv) an d which are again degenerate. Linearly combining these two 

states to two eigenstates of the rotation O, one obtains 




(4.59) 



Operating on the states ^±(r) with the 90 degrees rotation O, one now obtains the 
eigenvalues ±1, which implies that represents an s-wave, while \fr?_(r) actually 

has d-wave symmetry. As a consequence of an interplay of the various symmetries, the 
two states are exactly degenerate. The corresponding probability density is illustrated 
in figure 4 (right panel). Interestingly, although the states have different symmetries, 
their probability densities are identical. 

Under the discrete symmetries D' t , O, and R, the excited states *&±(r), which are 
bound by magnon exchange, transform as 




(4.60) 



States with d-wave symmetry can also be constructed for two holes localized on 

a Skyrmion. For example, for n = 1, the states _j i 2 o( r a, Xa, xp, 7) and 

^±0 2-11 o( r a, Xa, r p, Xp, 7) have d-wave symmetry. Under the symmetries D[ and D' 2 

they transform into #^Li ;1 - 2 ,o,o( r «' Xa, r p, Xp, l) an d *S-a,o ,-i,i,o( r «> Xa, rp, Xp, 7), 
which have s-wave symmetry. According to eq. (14.431) . under the various symmetries 
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the d-wave states transform as 

Di *±%,l,0,2,o( r ^Xa,r>X/3,7) = *^_i >ll _2 I 0,o( r a,Xa > ^,X^T). 

^VM,0,2,o( r «> Xa, ?> X/3, 7) = -*S-l,l,-2,0,o(»"a, Xa, ?> Xp, 7), 

°*?-l,l,0,2,o( r a,Xa,r /3 ,X/3,7) = -^S-l,l,0,2,o( r a>Xc" r /3 ) X/3 ) 7) ) 

H ^?-l,l,0,2,o( r « ; Xa,^,X/3,7) = *^0,2,-l,l,o( r a,Xa,?>X/?,7), 

2 ^*?0,2,-l,l,o( r a J Xa,r' i 8,X J 9,7) = *S-2,0,-l,l,o( r «i Xa, r /3, Xp, l) , 

fl ^S,2,-l,l,o( r m Xa, ^/3, X/3, 7) = -*^-2,0,-l,l,o( r a» Xa, ?> 7), 

°*±A2,-l,l,o( r a,Xa,r /3 ,X/3,7) = ~*± fo,2,-l,l,o( r a, Xa, ^ , X/3 , 7) , 

fl *S,2,-l,l,o('"a,Xa > r- i 9,Xi9 > 7) = *^-l,l,0,2,o( r a, Xa, Xp, 7)- (4-61) 

Similarly, the s-wave states transform as follows 

0i *±%,l-2,0,o( r «'Xa,^,X^,7) = *^_i,i,o,2,o( r «,X« s 77J,X0,7), 
2 ^*S-l,l,-2,0,o( r «> Xa, Tp, X/3, 7) = -*S-l,l,0,2,o( r a, Xa, Tfa Xp, 7), 

°*?-i,i,-2,o,o( r «'Xa,r /3 ,X/3,7) = *±!-i,i,-2,ofl{ra,Xa,rp,Xfi,i), 

R ^±^l,l,-2,0,o( r «'Xa,r /3 ,X/3,7) = *^-2,0,-l,l,o( r a,Xa,r /3 ,X/3,7), 

0/l *±%A-l,l,o( r «'Xa,^,X^,7) = *TA2-l,l,o( r a,Xa,'" / 8,X^,7), 

^*S-2,0,-l,l,o( r a, Xa, ?> X/3, 7) = -*S,2,-l,l,o( r «> Xa, ?> X/3, 7), 

°^?-2,0,-l,l,o( r a,Xa,r /3 ,X/3,7) = *±f-2,0,-l,l,0 ( r «> Xa, ^ , Xp , 7) , 

R ^± -2,0 -l,l,o{ r a,Xcc,rp,Xp,l) = ^-l,l,-2,0fl(Ta,Xa,rp,Xp,'r)- ( 4 - 62 ) 

Hence, just as for the magnon- mediated excited states, as a consequence of the inter- 
play of the various symmetries, s- and d-wave states are again degenerate. 

Alternatively, d-wave states also arise for two holes localized on a Skyrmion with 
winding number n = 2. For example, the two states ^± _i 2 -1 2 o( r a, Xa, fp, xp, l) have 
d-wave symmetry, and they transform into the states ^°Z-i\-2 i,o( r a, Xa, r p, Xp, 7), 
which again have s-wave symmetry, under D[ and D' 2 . According to eq. fl4.43p . the 
d-wave states transform as 

Di *if-l,2-l,2,o( r a, Xa, Tp, Xp, l) = *^-2,l -2,1,0 ( r a, Xa, Tp, Xp, l), 

D '*if-1,2 -l,2,o( r a, Xa, 7*0, X/3, 7) = -*5v-2,l,-2,l,o( r a> Xa, 7^, X/3, 7), 

^-l,2,-l,2,o( r a, Xa, 7-0, X/3, 7) = -*if-l,2,-l,2,o( r a, Xa, T"/?, X/3, 7), 

R y±,-l,2-l,2,o( r 'x>Xa,rp,Xp,l) = *J-l,2-l,2,o( r a,Xa,^,X/3,7), ( 4 -63) 
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while the s-wave states transform as 

D '*?-2,l -2,l,o( r a> Xa, rp, Xfi, 7) = *^-l,2,-l,2,o( r a> Xa> r /3, X/3, T)) 

2 ^*^-2,l,-2,l,o( r « J Xa, X/3, 7) = -*^-l,2,-l,2,o( r a, Xa, Tp, Xf3, 7), 

O *S-2,l,-2,l,0( r a,Xa,^,X/3,7) = ^±f-2,l,-2,l,0 ( r a , Xa , Tp, X/3, 7) , 

R ^_2,l,-2,l,o( r a,Xa,r/3,X/3,7) = ^-2,1,-2,1,0 (r a , Xa, Tp, \fa l) ■ (4-64) 

Depending on the details of the dynamics, whose investigation goes beyond the scope 
of the present paper, it may be possible that the degenerate s- and d-wave states have 
a lower energy than the p-wave states discussed earlier. 



4.7 Possible Implications for Cooper Pair Formation in High- 
Temperature Superconductors 

As we have seen, two holes, one of flavor a and one of flavor (3, can both get localized 
in the ground state of a single rotating Skyrmion with n — 1, which turns out to 
have p-wave symmetry. Alternatively, the holes may get localized on a rotating n = 1 
or 2 Skyrmion with s-wave or d-wave symmetry. As discussed in Appendix A, two 
holes of the same flavor can also get localized on a Skyrmion. It will be the subject 
of a subsequent publication to decide which of the various states is energetically most 
favorable. 

While in this paper we have concentrated on a detailed symmetry analysis, we also 
want to get at least a crude estimate of the binding energy of two-hole states localized 
on an n = 1 Skyrmion. Ignoring contact interactions between the two holes, the total 
energy of the bound state of two holes and a Skyrmion can then be estimated as 

E tot = 2M + 47rp s + 2E , (4.65) 

while two free holes (not localized on a Skyrmion) just have their rest energy 2M. 
Using the result of eq. fl4.13p . the perturbative vacuum thus becomes unstable against 
the formation of two-hole-Skyrmion bound states when 

4vrp s + 2E < 0.270MA 2 > Aixp s . (4.66) 

Hence, for sufficiently small spin stiffness p s , the instability will indeed arise. Similar 
instabilities are related to the formation of spiral phases in the staggered magnetiza- 



tion order parameter. In particular, in 25] we have shown that the ground state with 



a spatially constant staggered magnetization becomes unstable against the formation 
of a 45 degrees spiral phase for M'A 2 > 47tp s . Since antiferromagnetism is weakened 
upon doping, p s is expected to eventually go to zero. Before this happens, pairs of 
holes will get localized on a Skyrmion. 
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In order to get at least a rough idea of the involved energy scales, let us estimate 
the values of the relevant low-energy parameters for realistic lightly doped quantum 



antiferromagnets. By comparison with 23|, |26|, |27|, where a generalized t-J model on 



a square lattice with spacing a was considered at J/t ~ 0.3, one obtains the rough 
estimate 

M' « -L « A « 2.5Ja. (4.67) 

It would be interesting and definitely feasible to extract these parameters with high 
precision from numerical simulations. In this way, in the Heisenberg model (i.e. the 
undoped t-J model) very accurate numerical results have been obtained for the spin 
stiffness, the spinwave velocity, and the staggered magnetization per lattice site [l|,|3,|3| 

p s = 0.18081(11) J, c = 1.6586(3) Ja, M s = 0.30743(l)/a 2 . (4.68) 

Hence, one obtains 0.270M'A 2 « 0.5J compared to 47rp s = 2.2721(1) J, which implies 
that two-hole-Skyrmion bound states are still far from being energetically favorable 
at zero doping. The exchange coupling of undoped La2Cu04 is J = 1540(60) K [l|. A 
high transition temperature of T c m 50 K, and hence T c m 0.03J, would thus require 
a two-hole-Skyrmion bound state energy of about 

Atx Ps + 2E = 4vrp s - 0.270M'A 2 « -0.03 J p s « 0.04J. (4.69) 

If doping reduces p s by a factor of about 4 or 5 (and assuming for simplicity that the 
other parameters remain unchanged), the estimated energy scales should indeed be 
of the right magnitude in order to make two-holes localized on a rotating Skyrmion a 
viable candidate for a preformed Cooper pair of a high-temperature superconductor. 
Using eq. fl4.12p . one can estimate the radius of the Skyrmion, which sets the scale for 
the size of the candidate Cooper pair, as p ~ 1/(0.271M'A) w 5a, which again seems 
reasonable. 

It may involve some wishful thinking to assume that the d-wave state of two 
holes localized on an n = 1 or 2 Skyrmion will not only turn out to be energetically 
favorable, but also ready to condense at sufficiently large doping. However, we think 
that it is worthwhile to take this possibility seriously. Deciding whether the radial 
dynamics favors these states as promising candidates for a preformed Cooper pair 
in the pseudo-gap phase is the natural next step. The question of condensation is 
another important issue. 



5 Conclusions 



We have performed a detailed study of the localization of holes on a Skyrmion in an 
square lattice antiferromagnet. When two holes get localized on the same Skyrmion, 
they form a bound state. Interestingly, in some cases, the quantum numbers of these 
topologically non-trivial bound states are the same as those of the topologically trivial 
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bound states resulting from one-magnon exchange between two holes. The ground 
state of two holes weakly bound by one-magnon exchange has p-wave symmetry and 
may evolve into a strongly bound state of two holes localized on an n = 1 Skyrmion 
at strong coupling. 

Magnon-mediated two-hole bound states which are excited in the angular motion 
have s- or d-wave symmetry. Remarkably, s- and d-wave states are degenerate due to 
an interplay of the various symmetries. Similarly, there are strongly bound states of 
two holes localized on an n = 1 or 2 Skyrmion which also have s- or d-wave symmetry, 
and are again degenerate. Which of these states is energetically most favorable will be 
an interesting subject for future studies. If a d-wave state turns out to be the ground 
state at sufficiently strong doping, two holes localized on a Skyrmion are a promising 
candidate for a preformed Cooper pair in the pseudo-gap regime. Interestingly, the 
effective theory provides detailed predictions for the anatomy of these objects. In 
particular, their angular structure follows unambiguously from our symmetry analysis, 
and is insensitive to the details of the radial dynamics. 

Understanding the dynamical mechanism responsible for high-temperature super- 
conductivity has proved to be one of the most challenging problems in theoretical 
physics. While hole pair localization on a rotating Skyrmion may ultimately turn out 
not to be the relevant mechanism, it seems rather promising. Beyond the symmetry 
analysis presented here, studying its dynamics in more detail is certainly worthwhile. 
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A Hole Pairs of the Same Flavor 



In this appendix, we consider a pair of holes in a square lattice antiferromagnet 
residing in the same hole pocket. First, we investigate two holes localized on a rotating 
Skyrmion, and then we compare the results with the corresponding two-hole magnon- 
mediated bound states. 



A.l Schrodinger Equation for a Pair of Holes of the Same 
Flavor Localized on a Rotating Skyrmion 



Let us consider bound states of two holes of the same flavor / localized on a rotating 
Skyrmion. In this case, as a consequence of the Pauli principle, the holes cannot 
occupy the same quantum state. We distinguish the holes by an unphysical label 1 or 
2. In order to satisfy the Pauli principle, the wave function must be anti-symmetric 
under the exchange of the two labels. 



The Hamiltonian for two holes of the same flavor / is then given by 

H = H 1 + H 2 + W< 
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with H± ± and H±±±± given in eq. fl4.3ip . 



(A.l) 
/ 

(A.2) 



Before anti-symmetrizing the wave function in the artificial labels 1 and 2, we 
ignore the Pauli principle, and make the following ansatz for an energy eigenstate of 
two holes (distinguished by the labels 1 and 2) 

*^., mL ,„., mi , m (ri,Xl,r 2 ,»,7) = 

n, r 2 ) exp (ia [m\xi + ™\x-i ~ C/f ] ) exp(za(m - 1)7) ^ 
o"0"/V ; ffi m^mL,m^mi,m,+-( r i' r 2) ex P ( ia [ m \xi + mLx^ ) exp^cmry) 
^crf^a,m\,ml,m 2 + , m l, m -+(r 1 ,r 2 )exp (ia [m}_Xi + m+X 2 ]) ex.p(iarrvy) 
\ ^a,ml,m 1 _,ml,ml, m ,— ( r i,r 2 )exp (ia [mlxi + m -X2 + a f ^])exp(ia(m + 1)7) J 

(A.3) 
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As before, this solves the Schrodinger equation only if m l _ — m\ = n + 1, i = 1, 2. In 
this case, m is again an integer. The resulting radial Schrodinger equation now takes 
the form 



Hr^a,m\,m 1 _,m 2 + ,m 2 _,raiTl^ r 2) ^a,m} + ,m l _ ,ra 2 + ,m 2 _ ,m^cs,m} + ,m l _ ,ra 2 + ,m 2 _ ,m Vl i r 2/i (-^--4) 



with 



VV,mi ,m}_ ^m 2 ^ ,m 2 _ ,m,H — vl 1 ^2 ) 
^a,m^,m}_,rn 2 t _,m 2 _,m, — hv 1> r 2) 
\ i>a,m 1 + ,m 1 _,m 2 + ,m 2 _,m,— ( r l, r 2) J 



The radial Hamiltonian is given by 

H r = Hi + H^ + H^, 



(A.5) 
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HI J 

The matrix elements of the fermionic part of the radial Hamiltonian are given by 



H r++ 



1 



2M> 



D 2 
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= H*=V2A 
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r 2n + p 2n 



(A.8) 
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while the rotational Skyrmion contributions are given by 

n 2 ( P 2n p 2 " 



r++++ 2P(p) P 2 V 27T r 2n + p 2n r |n + p 2n 

n 2 ( p 2n p 2n x2 



Hr+-+- = , n rn + em- 7j—. — =- + 



r+ ~ + ~ 2V(p)p 2 V 2vr r 2n + p 2n r 2n + p 2n ' 

n 2 ( p 2n p 2n 

ii J. i i_ = . . - m + an- — h 



r ~ + ~ + 2V(p)p 2 V 2vr r 2n + p 2n r 2n + p 2n 

n 2 ( p 2n p 2n 

HI = — - - m + an h 1 + sr- — =- + — =- ] . (A. 9) 

r 2V(p)p 2 \ 27T r fn + p 2n r 2n + p 2n I 



A. 2 Symmetry Properties of a Pair of Holes with the Same 
Flavor Localized on a Skyrmion 



The spin operator / is again given by eq. (14.391) . such that 

,0i,Xi,r2,X2,7) = 



a,m\ ,mt_ ,m?^_ ,m 2 _ ,m v 

<^ ^ T ff 

m + an — W 
2tt J 



(ri,Xi ) r 2 ,X2 ) 7)- 



(A.10) 



Since m is an integer, at least for = 0, the state with two holes of the same flavor 
localized on a Skyrmion again has integer spin. 

The symmetries D' { , O, and R act on the two- hole wave function 

/ ^^ / „ i++ (ri,Xi,r 2 ,X2,7) \ 

1 K,n,-+{ri,Xi,r 2 ,X2,V 

V *£!»,— ( r l>Xl ; r 2,X2,7) / 



(A.11) 



as follows 



/ * 



D ^ai( r i> Xi> r 2, X2, 7) = exp(2ifc/a) 



(^l,Xl,^2,X2,7) \ 



*in 1 -+(n ) Xi ) r 2 ,x 2 ,7) 



vf^ / ni+ _(ri,Xi,r 2 ,X2,7) 
// 



-( r i ; Xi,^2,X2,7) / 



/ ^ / „,++(ri,Xi + f^2,X2 + f,7 



o 



*a / n( r l ) Xl,^2,X2,7) 



(T/^ n + _(ri, xi + f , r 2 , X2 + f , 7 - n%) 



a f ^ f J^_ + (r l} xi + f , r 2 , X2 + f , 7 



-( r i> Xi + f ,7-2, X2 + f , 7 - nf ) 



l ^ f Jn( r uXur 2 ,X2,l) 



2 5 ' Zi A.Z 1 2 

. ,++v i,-Xi,r 2 ,-X2, -7) \ 

*£fn,+-( r l, ~Xl> r 2, -X2, -7) 

*JJ lj _+(ri, -xi, r 2 , -X2, -7) 
V tfJJ*,— (ri,-xi,r 2 ,-X2,-7) / 



/ 



(A.12) 
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For the two-hole energy eigenstates this implies 

(n,Xi,r 2 ,X 2 ,7) ) 
2.,m 2 _,m( r i'Xi,^2,X2,7) = -exp (ier[mi + mi - mn) |) 

(n,xi,^2,X2,7), 



2 



a/m\ ,m 1 ,mP, ,m 2 ,m 
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,( r i ; Xi,^2,X2,7) =exp \ ia[m\ + ml - mnj- 

x*^i „i _ 2 m2 _(ri,Xi,r'2,X2 J T)> 



^ff,m|,m|.,m5.,m?.,m( ri ' Xl) r 2) X2; 7) ^ 



<r,m^j_ ,m_ ,mi_ ,m_ ,m \ 



ra\ .m 1 ,m 2 ,m 2 ,m 



(ri,Xi,^2,X2,7), 
(n,Xi,^2,X2,7)- (A.13) 



Here we have again assumed an appropriate phase convention for the radial wave 
function Vv.m 1 .m* , m \,m 2 ,m( r i) r 2)- in the context of the shift symmetries D[ we have 
used 



^.1 



m\ ,mL j*" 2 - ,™< 2 _ 
m 1 ^ ,m?_ jm^jjn^ ,m,- 



.(n,r 2 ) = 



112 2 



(ri,r 2 ; 



('"lj ^2) ^—a,—m}_ ,—m]^,—in 2 _ ,— m 2 ^ ,— m,H — vl ) ^2) > 



_ ,m 2 ^,m : L ,m,++ vl ' ^2) V' '— <r,— mL ,- 



.(r x ,r 2 ). (A.14) 



These relations follow from the symmetries of the radial Schrodinger equation (1A.4j) . 
In the context of the reflection symmetry R we have used 

^Pa,ml^ ,m\_ ,m^_ ,m?_ ,m,H — (- vl j ^2) ^P—a,ml^ ,m]_ ,m\ ,m?_ ,m,H — (- vl ? ^2) 5 

0CT,roi ! w|_ ,m3_ jnj 2 ^ ,?7i,H — ('"lj^j ^-(j,m| ,m?_ ,m^,m 2 ,m,H — ylt^2)t 

^a,m},,rn?_,mZ_,mL,m, — (-v l'^) '0— <r,mi_ ,mL 1WI4. iW 2 .. ,m, — (-vl> 

VV.mi. ,mL iTOj. ,m 2 ,m, v*l : ^2) '^P—a,m 1 ,,rn i _,m 2 i,m' 2 ^,m, vlj^)- (A. 15) 



The relations in eq. (1A. 15h follow from the symmetries of the radial Schrodinger equa- 
tion (IA.40 for = 0. As before, for G 7^ or tt, the Hopf term explicitly breaks the 
reflection symmetry. 

Let us now impose the Pauli principle by explicitly anti-symmetrizing the wave 
function in the artificial indices 1 and 2. For this purpose we act with the pair 
permutation P, i.e. 



^„_ + (r 2) X2 ) ri,Xi,7) 
.( r 2,X2,ri,xi,7) 



// 

OVf. 



(A.16) 



\ ^afn,— ( r "2,X2,?"i,Xl,7) / 
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For an energy eigenstate this implies 

F HU.^i,.iJ r - »■ *• »■ = *■ » »■ 7) - 

Here we have assumed a symmetric radial wave function, i.e. 



VV,?7ii_ ,m]_ ,m^_,m^ ,m,H — v 2. '"l) 

-( r 2,r"i) 



,m?_ ,m^_ ,mi_ ,?7i,H — h V 



2 2 11 



r ,m5_ ,m_ ,mi ,m_ ,m,+ 
"0(T,m5_ .m 2 ^ jm 1 , ,m, V 1) ^2)- 



<T,mi ,ml ,m^_ ,m?_ ,m 



(ri,Xi,r 2 ,X2,T) = 

(^l,Xl,^2,X2,7) 



(A.17) 



(A.18) 



(A.19) 



The properly anti-symmetrized wave function now takes the form 

^a / n( r l ) Xl,^2,X2,7) = -y= Ol . Xl , T 2 , X2 , l) ~ ? ^ S Jj?\ > Xl . r 2, X2 , 7)] 

For an energy eigenstate this implies 

Mr 112 2 

I 

As expected, in order to obtain a non-vanishing wave function, the two sets of quan- 
tum numbers mi, and m 2 + , m 2 .. must be different, because otherwise two identical 
fermions would occupy the same single particle state. If one would consider an anti- 
symmetric radial wave function, one could allow m\ = m\ and m\_ = m?_. 



(A.20) 



Based on eq. flA.13j) . the properly anti-symmetrized two- hole energy eigenstates 
transform as follows 



,m 2 _ ,m( ri ' Xl. r2 ' ^) ^-<r,-TO:L,--mi 



-m_ -ml-m 



(n,Xi,^2,X2,7), 



119 9 



7T 



o 



(n,Xi,^2,X2,7) = -exp ^a[m^ + - mn\- 

x *J^,mL,m» y_,m( r i> Xi, r 2 , X2, 7), 
2 1 



mn\ 



X 1122 



RiT,aa 



^ CT,m]j_,m^,m+,m?_,mV 1> Xl, r 2. Xli 7) ^ 



AS 

CT,m J ; /m", ,m" ,m 



Xri,Xi,r 2 ,X2,l)=^Z. 



119 2 



(n,xi,^2,X2,7), 
(ri,Xi,r 2 ,x 2 ,7), 

(n,Xi,^2,X2,7)- (A.21) 



In order to show this for the rotation O, we have used m} + + m 2 _ = m 2 + + m}_ 



Finally, let us combine states with flavors aa and (3(3 to eigenstates of O, 



± 



i7,m\ ,m?_ ,m?_ ,m?_ ,m 



(ri,Xi,r 2 ,X2,7) 
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71 



*r m i,m 1 _,m^mi,m( r l. Xl. r 2, X2, 7) ± ^J^.Tnl.m^mi.J 7 "!' Xl, r 2, *2, 7) , (A.22) 
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which transform as 

Di$ Jm|,mL, m 2 + y_,m( ri ' Xl > r2 ' X2 ' 7 ^ = ~^-a-ml -m\ ,-m 2 . ,-m% -m^ 1 ' Xl, r 2, X2, 7), 

°Kml,ml,ml, m l,J r ^ Xi, ^2, X2, l) = ±i exp (ia[m\ + m 2 _ - mn]|) 

X *Jml ,mL y + ,ml,m( r l' Xl, ^2, *2, 7), 
.raL.m^.ml.mC 7 "!' Xl, ^2, *2, 7) = .mi^.m^m^' Xl, r 2 , *2, 7)- (A.23) 

The lowest energy states in the same flavor channel are expected to correspond to 
m]_ = — l,m}_ = l,m\ = —2,m 2 _ = 0,m = or m} + = —l,m]_ = l,m+ = 0,m 2 _ = 
2, m = 0. These states transform as 

D **J-i,i - 2 ,o,o( r i> Xi, r 2 , X2, 7) = -*t ff ,-i,i, ,2,o( r i> Xi, r 2 , X2, 7), 

-2,0,0 ( r i> Xi, ?"2, X2, 7) = ±0-*J-i,i,-2,o,o( r i> Xi, r 2 , X2, 7), 
^,-1,1,-2,0,0 (^l, Xi, r 2 , X2, 7) = ±^^,-1,1,-2,0,0 Xi, r 2 , X2, 7), 

D ^J-l,l,0,2,o(n, Xl, r 2 , X2, 7) = -^a,-l,l,-2,0,o( r l' Xl, T 2 , X2, l) , 

°^-i,i,o,2,o( r i ; Xi, r 2 , X2, 7) = ^0-^^-1,1,0,2,0(^1, Xi, r 2 , X2, 7), 
^-1,1,0,2,0 ( r i, Xi, r 2 , X2, 7) = ±^-,7,-1,1,0,2,0 ( r i> Xi, ^2, X2, 7)- (A.24) 

This implies that the states *+ i _i j i,_ 2 ,o,o( r i, Xi, r 2 , X2, 7), *~ -1,1,-2,0,0 ( r u Xi, r 2, X2, 7), 
®- ,-i,i,o,2,o( r i,Xi,'"2,X2,7), 5+-i,i,o,2,o( r i,Xi,^2,X2,7) axe s-waves, while the states 

J+ -1,1,0,2,0 ( r l, Xl, ^2, X2, 7), *I,-l,l,0,2,o( r l, Xl, ^2, X2, 7), *i,-l,l,-2,0,o( r l, Xl, ^2, X2, 7), 

*+,-i,i,-2,o,o( r i,Xi,^2,X2,7) are d-waves. 

A. 3 Comparison with Magnon-Mediated Two-Hole Bound 
States of the Same Flavor 



In [14| states of two holes of the same flavor bound by one-magnon exchange have also 
been investigated. Here we summarize as well as extend some of the relevant results. 
We consider two holes of the same flavor / with opposite spins + and — . In the rest 
frame the wave function depends on the distance vector f which points from the spin 
+ hole to the spin — hole. Since magnon exchange is accompanied by a spin-flip, 
the vector f changes its direction in the magnon exchange process. The Schrodinger 
equation thus takes the form 

- -La# (r) + V f, (r)V(-r) = E^(f). (A.25) 

The one-magnon exchange potential for two holes of the same flavor is given by 

v A^MM ^ (0 = __*L*^. (A.26) 

2np s r 2 2np s r 2 
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Figure 5: Angular wave functions ceo(ip — j, \M'^) (solid curve) and sei((p — j, \M'^) 
(dashed curve) as well as angle- dependence sin(2<£>) of the potential (dotted curve) for 
two holes of flavor a residing in a circular hole pocket (M'A 2 /4Tcp s = 1.25). 

We make a separation ansatz 

*(f)=#(r)x'(p). (A.27) 
The ground state is even with respect to the reflection of r to — r, i.e. 

X n (<p + n) = x n (<p)- (A.28) 
The angular part of the Schrodinger equation then takes the form 

- ± ^ M*P)7&{<P) = -Aix'i(^)- (A.29) 

Here + and — are associated with an aa and a (3(3 pair, respectively. Again, eq. (1A.29j) 
is a Mathieu equation. The ground state with eigenvalue — Ai takes the form 

n . , -i . 7T. 1 . 7T M'A 2 , 

X ± M = xA<P - 3 ) = ^ce„(^ - j.i^), 
A1= UP!V + 0(A8) . (A . 30 ) 



2 V 4vrp s 

The first excited states are odd with respect to the reflection of r to — r, i.e. 

x2(v + ir) = -x2(*0, (A.3i; 
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and the angular part of the Schrodinger equation now reads 

T sm(2<p) X °(<p) = -Aax'J^)- (A.32) 



d 2 X 'Uv) M'A 2 



Now — and + are associated with an aa and a (3(3 pair, respectively. The excited 
states with eigenvalue — A 2 are given by 

a, \ 2/ *\ 1 . 7T M'A 2 . 



,,, ,, ,i i 7T M'A 2 x 





l 




=ce 








'M'A 2 \ 







M'A 2 1 M'A 2 \ 2 1 /M'A 2 \ 3 ^ /AS , /A . 
A 2 = -l + + - +CA 8 . A.33 

The angular wave functions for the ground state and for the first excited state together 
with the angular dependence of the one-magnon exchange potential are shown in figure 
5. 

As before, the radial Schrodinger equation takes the form of eq. (14.531) . Again, the 
short-distance repulsion between two holes is modeled by a hard core of radius r' , i.e. 
R'{r' Q ) = 0. The value of r' Q may, however, differ from r in the a(3 case. The radial 
wave functions are thus given by 



R' i {r)=A' t K u UM'\E' k \r) 1 = 1, 2, 3, . . . , u = iy/\, (A.34) 



and the energy is determined from K V {^J M'\E' ik \r' ) = 0. 

There are two degenerate states — one for an aa and one for a (3(3 pair, which 
are eigenstates of flavor related to each other by a 90 degrees rotation. The two 
degenerate states can be combined to eigenstates of the rotation symmetry O. For 
this purpose, we construct the 2-component wave functions 

= K(r) ( ) , *J« = «i(r) ( ) , (A.35) 

whose first component represents the aa and whose second component represents the 
(3(3 pair. Under the various symmetries, the two degenerate ground states transform 
as 



= JH(r) ( ) = fli(r) ( ± ^ ) = ±»3(0. (A.36) 



44 



Figure 6: Probability distribution for bound states of two holes with flavors aa or (3/3, 
combined to an eigenstate of the 90 degrees rotation symmetry O. Left panel: the 
ground state with p-wave symmetry. Right panel: excited states with s- or d-wave 
symmetry, but with identical probability densities (M' K 2 / A~Kp s = 1.25, r' = a). 

Again, the corresponding eigenvalues of the 90 degrees rotation O are o = ±i, and 
hence, as for a(3 pairs, the symmetry is actually p-wave. Similarly, the two degenerate 
first excited states transform as 



Again, the first excited states transform as s- or d-waves. The resulting probability 
distributions, which resemble d xy symmetry, are illustrated in figure 6 for the ground 
state (left panel) and the first excited state (right panel). Unlike for an a (3 pair, in 
the same flavor case the lowest energy bound states localized on a Skyrmion have a 
different transformation behavior than the magnon-mediated two-hole bound states. 




( Tx' 2 &) 



) 



±*'|(r). (A.37) 
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